EPJ manuscript No. 

(will be inserted by the editor) 



Minimal conductivity in bilayer graphene 

M. I. Katsnelson 

Institute for Molecules and Materials, Radboud University Nijmegen, 6525 ED Nijmegen, The Netherlands 
Received: date / Revised version: date 

Abstract. Using the Landauer formula approach, it is proven that minimal conductivity of order of e} jh 
found experimentally in bilayer graphene is its intrinsic property. For the ceise of ideal crystals, the conduc- 
tivity turns our to be equal to per valley per spin. A zero-temperature shot noise in bilayer graphene 
is considered and the Fano factor is calculated. Its value 1 — 2/7r is close to the value 1/3 found earlier for 
the single-layer graphene. 
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It has been observed recently that bilayer graphene, 
that is a two-dimensional allotrope of carbon formed by 
two graphite atomic sheets, has a minimal conductivity 
of order of jh T. The same property has been found 
earlier in the single-layer graphene . Both single- and 
bilayer graphene are gapless semiconductors, with conical 
and parabolic touching of electron and hole bands, re- 
spectively PEIEI- The charge carriers in the single-layer 
graphene are massless Dirac fermions which is a crucial 
point when explaining the conductivity minimum [311321 ■ 
Actually, this anomalous property of the two-dimensional 
massless fermions was considered theoretically be- 
fore discovery of graphene. A crucial physical phenomenon 
here is the Zitterbewegung of quantum ultrarelativistic 
particles 4 which plays a role of "intrinsic" disorder; the 
latter being confirmed by calculations of the shot noise in 
ideal graphene for zero doping which turns out to have the 
same value (Fano factor 1/3) as disordered metals jS]- At 
the same time, observation of the finite minimal conduc- 
tivity in bilayer graphene is a serious challenge for theory 
. Here I present a solution of this problem based on the 
same Landauer formula approach which was used earlier 
for the single-layer case j^IS]. 

The bilayer graphene is a zero-gap semiconductor with 
parabolic touching of the electron and hole bands described 
by the single-particle Hamiltonian [Tllll)| 

H=( ^ -{p,-^pyf/2m\ 
\ - [px + WyT I'^m / 

where pi — —ihd/dxi are electron momenta operators and 
m is the effective mass (here we ignore some complica- 
tions due to large-scale hopping processes which are im- 
portant for a very narrow range of the Fermi energies 10 ). 
Two components of the wave function are originated from 
crystallographic structure of graphite sheets with two car- 
bon atoms in the sheet per elementary cell. There are two 
touching points per Brillouin zone, K and K'. For ideal 




Fig. 1. Geometry of the sample. 



crystals, no Umklapp processes between these points are 
allowed and thus they can be considered independently. 
Our final result for the conductivity should be just multi- 
plied by four due to two touching points and two spin pro- 
jections (we will not take into account electron spin explic- 
itly in our consideration) . To calculate the conductivity at 
zero energy we will use the Landauer formula expressing 
the conductance of the system in terms of transmission 
coefficients. Similar to Ref. @] we will use the simplest 
boundary conditions assuming that the sample is a ring 
of length Ly in the y-direction and leads connected with 
the sample at x = and x — are made from doped 
bilayer graphene with potential Vo > and Fermi energy 
Ef = ~Vo = -h^k'j,/2m (Figure 1). 
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Let us first find tlie solution of tlic Sclirodinger equa- 
tion with zero energy, = where is a "spinor" with 
components tpi and ip2- They satisfy the equations 



d 9 , , „ 

^ * «~ ■02 = 0, 

ax oy/ 
ox oy/ 



(2) 
(3) 



Due to the periodicity in the y direction both wave func- 
tions should be proportional to exp (ikyy) where ky — 



2Tm/Ly 



0,±1,±2,.... This immediately gives us the 



following x-dependence for the wave functions: 

V-i (x) = (^lX-HBl)e'=«^ 
V'2 (x) = {A2X + B2) e-^y^ 



(4) 



{0 < X < Lx). The constants Ai and Bi should be found 
from the boundary conditions at a; = and x = which 
are nothing but continuity conditions for both functions 
tpi and V'2 and their derivatives [TT] . 



It will be shown further that the values of ky essential 
for the electron transmission are of the order of Ly^ and 
thus much smaller than kp in the leads. Therefore, one 
can restrict ourselves to the case of normal incidence only 
for the wave functions outside the sample: 



ipi (x) = e'^"'' +re-"^' -^ +ce 



-ikpx 
-ikpx 



kpx 
kpx 



(5) 



for X < and 



ipi {x) = te''=^(^"-^=) + de-'"'^''-^''\ 
1P2 {x) = te^'^'^i''-^-) _ de-'"''^''-^''^ 



(6) 



one can easily solve the equations lO and find 
2iLx cosh{kyLx) 



for x > Lx- Here r and t are reflection and transmission 
coefficients, respectively. One should stress that to satisfy 
all the boundary conditions for the case of a bilayer, one 
has to include not only oscillatory but also exponentially 
decaying solutions of the Schrodinger equation [111. 

Using the boundary conditions at the sample-lead bound- 
ary, one finds the set of linear equations 



B2, 

Al+Blky, 

A2 — B2ky, 
t + d, 
t-d, 
k {it — d) , 
{-F2ky + A2)X-' = k{it + d) 



i + r + c ■ 
1 + r — c 
kp [i (1 - r) + c] : 
kp [i (1 — r) — c] : 
FiX ■ 

A,)X. 



(Fiky 



where F.^ = AiL^ + Bi,X = exp (kyLx) 
By use of the assumptions 

kp ^ ky , F,^ ^ 



(7) 



(8) 



t = 



kp Ll + ^cosh^{kyLx)' 



(9) 



Corrections to this formula are of order of l/{kpLx); we 
cannot keep them in the answer since terms of the same 
order of magnitude have been omitted by considering the 
normal-incidence case for the wave functions in the leads 
®, ®. 

Thus, for the transmission coefficient T„ = \t (ky — iim/ L 
one obtains the final result 



ikjpLl cosh^ (kyLx) 
k%Li + Acosh^ [kyFx) 



(10) 



One can see that the transmission coefficient reaches 
the maximum value equal to 1 at cosh (kyF^) ~ kpF^/V^ 
or, approximately, at 



ky/kp ~ In ( V2kpLx ) / (kpLx) 



(11) 



which obviously satisfies the condition ^ for macroscop- 
ically large F.^ ^ k^^. Note that the complete transmis- 
sion through the potential barrier for some finite incident 
angles is a characteristic property of the bilayer case, in 
contrast with the single layer, where the complete trans- 
mission takes place at exactly normal incidence |11| . 

Using the Landauer formula (for review, see Refs. |12l 
I13p one can calculate the conductance per valley per spin 



(12) 



n— — 00 



Similar to Refs. to calculate the conductivity of bi- 
layer graphene at zero energy one should consider the case 
Fy ^ Fx- In that case the sum in Eg. 1)12(1 can be re- 
placed by an integral. Introducing the integration variable 
z = cosh {2kyFx) + l and taking into account the condition 
(|HJ one finds for the conductivity a = {F^/Fy) g: 



2h' 



(13) 



Thus, the conductivity of bilayer graphene has the same 
order of magnitude than for the single-layer case (where 
the coefficient I/tt, instead of 1/2, was obtained by similar 
method in Refs. 0|S]). This result looks rather unexpected 
since the electron spectra in these two cases are drastically 
different. More accurate calculations of the integral gives 
a correcting multiplier 1 -I- + ... in Ea. (|13|l . 

Following Ref. |5j one can estimate the Fano factor 
characterizing the intensity of electron shot noise: 



F : 



E r„(i-T„) 

,— — 00 

CSO 



(14) 
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(for a general review of the quantum-limited shot noise 
and physical meaning of the Fano factor, see Refs. [I'M 
I14p. A straightforward calculation for the case Ly ^ 
gives us the answer 

F = l-- (15) 

TT 

which is rather close to the value 1/3 found for the case 
of the single layer, as well as for the case of disordered 
metals [H]- This means that, in a sense, the case of bilayer 
graphene is also characterized by some "intrinsic" disorder 
similar to the Zitterbewegung ^j. 

Unfortunately, the accuracy of experimental data ^ is 
not sufficient to establish the numerical coefficient in the 
expression for the minimal conductivity. For the case of a 
single layer it is close to 1, instead of I/tt. 

To conclude, we have demonstrated that Dirac energy 
spectrum is actually not important for existence of min- 
imal conductivity in graphene. The latter has the same 
order of magnitude both for conical (a single-layer case) 
and for parabolic (a bilayer case) energy spectrum near 
the band crossing point. 
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